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The role of the interaction with the nearest electronic state 3 Sj"_ on the hyperfine structure and 
magnetic properties of the a(l)[ 3 E^] state of PbO is assessed. The accounting for this contribution 
leads to difference between g- factors of the J = 1 il-doublet levels, Ag — 37 x 10~ 4 , that is in 
a good agreement with the experimental datum Ag — 30(8) x f0 -4 . The contribution of this 
interaction rapidly grows with J. For J = 30 the difference of ^-factors of f2-doublet states reaches 
f00%; for hyperfine constants it is 18%. These differences also depend on the electric field and 
for E — 11 V/cm for 207 PbO the difference in g-factors turn to zero. The latter is important for 
suppressing systematic effects in the electron electric dipole moment search experiment. 



The use of a(l) excited state of PbO molecule has been 
proposed to search for electric dipole moment (EDM) of 
the electron This experiment is a serious test of the 
"new physics" beyond the Standard Model including dif- 
ferent supersymmetric models • Because of that the 
molecule was intensively investigated both theoretically 
and experimentally [lll - 



In the adiabatic approximation rotational levels of the 
a(l) state of PbO are determined by the effective spin- 
rotational Hamiltonian 



H sr = B'J' 2 + Ay (J e • n)(I • n) + 
Ai B G||(J e -n)(B-n)-Z?n-E 



(1) 



Here B 1 is the rotational constant, J, J e , I are the 
electron-rotational, electron and nuclear angular momen- 
tum operators, respectively (in this paper we will mea- 
sure angular momentum in units of h) , E and B are exter- 
nal electric and magnetic fields, D is the molecular-frame 
dipole moment, n is a unit vector along the molecular 
axis, £, directed from Pb to O, /^b is Bohr magneton. 
The hyperfine constant A» and g-factor Gn are deter- 
mined by the expressions [14j 
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where S e is the electron spin operator, /ipb is the mag- 
netic moment of 207 Pb, a.i are the Dirac matrices for the 
i-th electron, is its radius- vector in the coordinate sys- 
tem centered on the Pb atom, £1 = (5 , 3 E + | J^|^ 32 + } = 

±1. From naturally abundant isotopes only 207 Pb (/ = 
1/2) has nonzero /i Pb , for 208 Pb and 206 Pb 1 = and, 
therefore, fipb = 0. 

The parameters B 1 , An, Gn and D can be obtained 
both theoretically from calculation of the electronic wave- 
function ^a(i) ano - fitting the experimentally observed 



transitions to the parameters of the spin-rotational 
Hamiltonian ([T]). Comparison of theoretical and experi- 
mental values gives us information about accuracy of the 
calculated wavefunction "^(l) ana -j therefore) also gives 
information about accuracy of the calculated effective 
electric field, Wd, seen by an unpaired electron 0, H|. 
Note, that Wd can not be measured independently, but 
it is required for extracting d e from the EDM experiment. 
The experimentally observed parameters A|| , Gn also can 
be used for a semiemperical evaluation of Wd [6] ■ 

Previous investigations of PbO were based on the spin- 
rotational Hamiltonian (JXJ) . The main goal of the present 
work is to account for the interaction with the nearest 



electronic state 3 E„ 



which modifies the form of this 



Hamiltonian. To the best of our knowledge this is the 
first investigation of such kind for open shell diatomics. 

In the present paper the hyperfine structure of rota- 
tional levels was obtained by numerical diagonalization 
of the Hamiltonian in the basis set of electronic rotational 
wavefunctions 



3V + 
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Oit,±da,0)U Ml ,*lrt- e J M Ja,/3)U Ml , (4) 



where J Mn {a,p) 



V(2J+l)/47r£>^(a,A7 = 0) 
and Umj are rotational and nuclear spin wavefunctions, 
M and Mj = ±1/2 are projections of the angular mo- 
menta, J and I, on the laboratory axis z. When elec- 
tronic matrix elements are known then matrix elements 
on the basis set (@|) can be calculated with the help of the 
angular momentum algebra [l5l ] . Required diagonal elec- 
tronic matrix elements, being, in fact, the parameters of 
the spin-rotational Hamiltonian ([1} , are known from ex- 
periments. For the fifth vibrational level of the a(l) state 
of PbO they are B' = 0.235296 car\A\\ = -4.1 GHz, 
G|| = 1.857, D = 1.28 a.u. \£Mm- For purposes 
of the present study, it is not required to account for 
the small difference between the rotational constants of 
206,207,208pj- ) o molecules. The differences in properties 
discussed below are relevant only to the fact that the 
isotope 207 Pb has hyperfine structure. The off-diagonal 
electronic matrix elements were calculated in the present 
study by the configuration interaction method with the 
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generalized relativistic effective core potential [la . Il7 |. 
The scheme of the calculation is the same as that in the 
paper @. The calculated matrix elements are 



^ = S'(*; E+ |^|*S E+ > =0.17 cm- 



(5) 



— (* e + l 
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= -0.7 GHz, (6) 



G± = 



3 e; 



I J! 



S£|#S E+ > = 1.5. 



(7) 



It is known that Hamiltonian (JTJ leads to two-fold de- 
generacy of levels with different signs of Q. This degen- 
eracy is in fact only approximate. When the interaction 
(O is taken into account each rotational level splits on 
two sublevels, called f2-doublet levels. One of them is 
even (p = 1) and the other one is odd (p = —1) with 
respect to changing the sign of electrons and nuclear co- 
ordinates. In order to reproduce experimental value of 
the ^-doubling, 5.6 J(J+ 1) MHz [12j], the matrix el- 
ement ((SJ has to be equal to 0.15 cm" 1 . We consider 
this a good agreement, but will use experimental value 
hereafter. The states with p = (— 1) J are denoted as e 
and with p = (— 1) ,/+1 as / states. Note that the wave- 
functions *3 S + J M (ct, {3)Um! are / states, and they do 

o- ' 

not interact (see below) with e states of the a(l), unless 
parity is not conserved, due to weak interactions. 

Interactions JHJl and j7]) lead to different hyperfine 
structure and magnetic properties of the e and / lev- 
els. One can estimate from the second order perturbation 
theory that contribution from the terms 



l(*3 S + #j\/,±lt / M / |-ffhfs(mag)| , I'? s + 
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is small. Here -ffhfs and -ff mag are Hamiltonians of the hy- 
pcrfine interaction and the interaction with the external 
magnetic field, respectively. However, the terms 

2Re((^ + ei L±1 U Ml \2B'33 e m^ 6 J M>Q U Ml ) x 

(*3 S + 6, M,O^M / |i?hfs(mag)|*3 S + 9 J M±1 U Ml )) 



are much larger and their influence on the spectrum of 
the a(l) state is observable. 

In TableUwe give calculated g-factors for / states with 
different quantum numbers J of 206 ' 208 PbO molecule. 
For e states calculated g e = 1.85700 and is independent 
on J. We define <?-factors so that the Zeeman splitting 



TABLE I: The ^-factors for /-states of 2 06,208 pbo ag a func _ 

tion of J. For e states g e = 1.85700 and is independent on 
J. 

J 9f J 9f 

1 1.86074 10 2.06255 

2 1.86822 12 2.14848 

3 1.87943 15 2.30537 

4 1.89438 20 2.64142 
6 1.93549 25 3.07055 
8 1.99155 30 3.59256 



TABLE II: Calculated values of HFS as a function of J for 
/- and e-levels of 207 PbO. Results of Ref.QJ] 

are given in the 

last column (see Eqs.(l) and (2) therein). These results are 
for both /- and e-levels. 

J this work Ref. [llj 
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is equal to g e ^f)psB z M / J(J + 1). Our calculations were 
done using the finite field method. 

The obtained difference gf~g e = 37 x 10~ 4 for J = 1 is 
in good agreement with the experimental result 5/ — g e — 
30(8) x 10" 4 [11]. As it is seen from TableQ]the difference 
is rapidly increasing with J, and for J = 30 gf is about 
two times larger than g e . Another point to note is that 
matrix elements (JS|) and ([7]) do not contribute to g e and 
it remains J-independent and unchanged. This is due to 
the mentioned above parity selection rule. Limiting by 
the terms © we obtain 



9f~9e 



A ■ Gj_ ■ J(J+ 1) 



E 



E: 



= 1.87032 x 1Q- 3 J{J + 1) 



that is in a good agreement with Table fl] 

In Table HU the hyperfine splitting (HFS) calculated 
between F = J — 1/2 and F = J + 1/2 levels as a func- 
tion of J is given for e and / states of 207 PbO. Also the 



results obtained by applying Eq. (1) and (2) of ref. 11] 
are listed. Eqs. (1) and (2) of ref. [1JJ give HFS in the 
framework of the Hamiltonian {]]). The interaction with 
the 3 T,Q- is not taken into account in the ([T]), therefore 



(9) Eqs. (1) and (2) of ref. 



ll| give the same HFS for e and 
/ states of the a(l). Similarly to g factors, the hyperfine 
structure of e states is not affected when interactions (O 
and ([6]) are taken into account. However, there is a small 
difference between the hyperfine splittings calculated by 
Eqs. (1) and (2) in Ref. [ll| and that calculated for e 
states in this paper. This difference is related with the 
fact that the mixing between the states with AJ = ±1 



3 



in Eqs. (1) and (2) of Ref. ll| is taken into account 
in the framework of the second order perturbation the- 
ory, whereas in the current work it is calculated more 
accurately, by using the numerical diagonalization of the 
Hamiltonian. 

In the electron EDM search experiment the Stark split- 
ting between J — 1,M ± 1 states of the e or / levels is 
measured. This Stark effect induced by the interaction 
with the electron EDM that violate both parity (P) and 
time reversal (T) invariance, and is not related with the 
(large) dipole moment D presented in the (HJ. For details 
see pp. 1-3 in Ref. j3|. In the external electric field the 
states J = 1, M = ±1 remain degenerate, unless both P 
and T are violated. However an external magnetic field 
remove degeneracy between them and can mimic the ex- 
istence of the EDM. For J = 1 levels the systematics due 
to spurious magnetic fields can be suppressed if the dif- 
ference between g e and g / can be made smaller [ijj • The 
external electric field mixes e and / levels. Therefore, on 
the first glance, one can expect that when increasing the 
electric field the initial small difference between q e and 
gf can be made zero. However, it was found in [12j that 
this difference for 206 . 208 pbO is actually increases as the 
electric field increases. This fact was explained by M.G. 
Kozlov (see acknowledgments in 13]) by accounting for 
the mixing with J = 2 level. In the present paper we 
reproduce this result for spinless isotopes of led and also 
calculate g-factors for J = 1, F = 1/2,3/2 states of the 
207 PbO. For 207 PbO, ^-factors was defined so that the 
Zeeman shift is given by 



F(F + 1) + J(J+l)-3/4 
2F(F + l)J{J + 1) 



Mp. 



With this definition they will coincide with <?-factors of 
206,208 PbO in the limit of zero hyperfine interaction. The 
corresponding results are given in Fig. (TIJ. One can see 
that difference between g e and gf for J = 1,F = 3/2 
does not converge to zero as E increases. However, for 
F = 1/2, J = 1 at E re 11 V/cm g e and gf become 
equal. The plotted g e and g f for J = 1 206,208p bo are 
in agreement with Fig. (5) of ref. [ijjj]. Large deviation 
of ^-factors for J = 1,F = 3/2 of 207 PbO from those 
for J = 1 of 206 > 208 pbO is explained by mixing of the 
J = l,F = 3/2 and J = 2,F = 3/2 levels of 207 PbO 
that is induced by the hyperfine interaction. 

In the EDM experiment the maximum Stark split- 
ting, 2W d • d e , between F = 1/2, Mp = ±1/2 lev- 
els is achieved for the fully polarized molecule. In 
Fig. ([2|) we plot the calculated Stark splitting between 
F = 1/2, Mp = ±1/2 levels as function of the external 
electric field. For E = 11 V/cm the obtained splitting is 
about 75% of the maximal value. 

In this work we account for non-adiabatic interaction 
of a(l)[ 3 £^~] state only with the state 3 E^|"_ . There are 
several reasons for this. One can see (lil . [l9| that the 
3 S^"_ state is the nearest one to the a(l) state. All other 



states, except 3 A, are more than an order of magnitude 
further away. Accounting for the non-adiabatic interac- 
tion with the 3 Ai state (the same = 1 as in a(l)) will 
lead only to a small modification of the parameters of the 
spin-rotational Hamiltonian (fT]). Since we use the exper- 
imental data, those interactions with the 3 Ai and other 
f2 = 1 states are taken into account. Though the interac- 
tion with 3 A2 can not be described in the framework of 
the Hamiltonian (fT]), it will not lead in the leading order 
to the difference in properties of the / and e states that 
is a topic of this paper. Moreover, our calculation show 
that the corresponding matrix element 

B'(^ t \J-m Af ) « 3 x 1(T 3 cm" 1 



3 AJ 



is small as compared to ([5]) . fl — 3 states are not mixed in 
the leading order due to the selection rule. The validity of 
the above approximation is approved by the fact that the 
calculated and the experimentally obtained differences 
of the g-factors for e and / J = 1 states are in good 
agreement. 

Finally we have investigated the influence of the in- 
teraction with the nearest electronic state 3 £,}- on 
the hyperfine structure and magnetic properties of the 
a(l)[ 3 £^~] state. We have shown that it is required for 
its accurate description, especially for g-factors. One can 
suppose that similar situation takes place also for other 
diatomics in CI = 1 states. It is found that the differ- 
ence between g e and gf for 207 PbO is converged to zero 
at E ks 11 V/cm. The latter is important for the sup- 
pressing systematic effects in the electron EDM search 
experiment. 
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FIG. 1: Calculated (/-factors for e(g e ) and /((//) states, 
(a) Solid lines correspond to J = 1,F = 1/2 hyperflne 
levels of 207 PbO, dashed lines correspond to J = 1 ro- 
tational levels of 206 ' 208 PbO. (b) Solid lines correspond to 
J = 1,F = 3/2, \Mf\ = 3/2, dashed lines correspond to 
J = 1,F = 3/2, \M F \ = 1/2 hyperfme levels of 207 PbO 
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FIG. 2: EDM induced Stark splitting between M F = ±1/2 
levels of J = 1, F = 1/2 state of 207 PbO 



